Cellular Dynamical Mean Field Approach to Strongly Correlated Systems 
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We propose a cellular version of dynamical-mean field theory which gives a natural generalization 
of its original single-site construction and is formulated in different sets of variables. We show how 
non-orthogonality of the tight-binding basis sets enters the problem and prove that the resulting 
equations lead to manifestly causal self energies. 



Dynamical mean field theory (DMFT) has been very 
successful in describing many aspects of strongly cor- 
related electron systems and presently much effort 
is put into implementing it for realistic calculations of 
materials properties of solids. By construction, this 
method describes correctly local correlations but misses 
altogether the effect of short-range order. This can be 
understood by noticing that in the limit of large lattice 
coordination (z — > oo) Q where the single-site DMFT is 
exact, the interactions which induce short-range correla- 
tions, such as magnetic superexchange, are scaled as ^ 
and in the absence of frozen correlations with z neighbors 
they disappear from the problem. The effects of magnetic 
correlations on single-particle properties are captured by 
single site DMFT in magnetically ordered phases HQ. 
However many applications require generalizations of the 
DMFT to capture short-range correlations, in the ab- 
sence of broken symmetries. This is an active area of 
research and several methods have already been put for- 
ward Jl],||-||] for this purpose. 

The intuitive idea behind cluster methods is to treat 
certain local degrees of freedom (cluster degrees of free- 
dom) exactly while replacing the remaining degrees of 
freedom by a bath of non interacting electrons which hy- 
bridize with the the cluster degrees of freedom so as to 
restore the translation invariance of the original prob- 
lem. The simplest example of this idea is the Bcthc 
Peierls cluster. The basic difficulty of a cluster tech- 
nique is determining a suitable self consistency condi- 
tion without generating unphysical solutions which vio- 
late causality. For a discussion of this longstanding prob- 
lem in the context of disordered systems see [Q . Ingersent 
and Schiller [Q and independently Georges and Kotliar 
[jlj , introduced a truncation of the skeleton expansion in 
real space, which can be represented by solving coupled 
impurity models of different sizes. This method is not 
manifestly causal, but recent work 0, suggests that the 
problems with causality encountered in the earlier treat- 
ments are the result of inaccurate approximations in the 
solutions of the impurity models. Jarrell || and collabo- 
rators have suggested an alternative cluster scheme, the 
dynamic cluster approximation (DCA) which is a cluster 
scheme in momentum space, whereby the cluster con- 
sidered, if regarded in real space has periodic boundary 



conditions. This body of work j6[ and its extensions |8|] 
established the computational feasibility and the exis- 
tence of an a priori causal cluster schemes. 

In this paper we pursue generalizations of the single 
site DMFT inspired by analogies with electronic struc- 
ture methods. This cellular DMFT (CDMFT), remains 
close in spirit to the DMFT ideas described in the intro- 
duction, where the clusters have free (and not periodic) 
boundary conditions. We prove two central points: 

I) The DMFT construction Q can be carried out in a 
large class of basis sets. This observation frees us from 
the need to introduce sharp boundaries in real space. 
This approach is inspired by ideas from electronic struc- 
ture, in which one achieves a cellular description by 
means of orbitals which can have a variable spatial ex- 
tension. 

II) This CDMFT construction is manifestly causal, 
i.e. the self energies that result from the solution of the 
cluster equations obey ImS(fe, uS) < 0, eliminating a pri- 
ori one of the main difficulties encountered earlier in de- 
vising practical cluster schemes. 

It is useful to to separate the three essential elements 
of a cluster scheme (see Fig. [3]): 

a) The definition of the cluster degrees of freedom, which 
are represented by impurity degrees of freedom in a bath 
described by a Weiss field matrix function Go- The so- 
lution of the cluster embedded in a medium results in a 
cluster Green's function matrix and a cluster self energy 
matrix. 

b) The expression of the Weiss field in terms of the 
Green's function or the self energy of the cluster, i.e. the 
self-consistency condition of the cluster scheme. 

c) The connection between the cluster self energy and the 
self energy of the lattice problem. The impurity solver 
estimates the local correlations of the cluster, while the 
lattice self energy is projected out using additional infor- 
mation, i.e. the periodicity of the original lattice. 

Our construction applies to very general models for 
which a lattice formulation naturally appears. It can be 
thought of as an extension of the band structure formal- 
ism that takes into account the electron electron interac- 
tions. The lattice Hamiltonian, H[fi a , fj a ], (one example 
could be the well known Hubbard Hamiltonian) is ex- 
pressed in terms of creation and annihilation operators 
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fiu and fj where i runs over the sites of a d dimen- 
sional infinite lattice i = (ii, . . . , id), the index a denotes 
an internal degree of freedom such as a spin index or 
a spin-orbital or band index if we consider an orbitally 
degenerate solid. 

a) Selection of cluster variables: The first step in a mean 
field approach to a physical problem, is a selection of a fi- 
nite set of relevant variables. This is done by splitting the 
original lattice into clusters of size 11^=1^ arranged on a 
superlattice with translation vectors R n . On this super- 
lattice we choose wave functions \R n a) partially localized 
around R n with a — 1, . . . , N denoting an internal clus- 
ter index. The relation between the new wave functions, 
\R n a), and the old ones, \io~), is encoded in a transforma- 
tion matrix, S Rn a^, SUCh that \R n a) = J2ia \ i(J ) ^ia,R n a- 

Due to the translation symmetry of the lattices we have 
S , j?„Q i i (T = S aer (r(i) — R n ) where r(i) is the position of 
site i. The creation and annihilation operators of the 
new basis are related to the operators of the old ba- 
sis by cr„ q = Ylicr ^Rna,i<rficr and the operators that 
contain the "local" information that we want to focus 
our attention on are c a = Cfu n =o)ai i- e - the operators 
of the cluster at the origin. We will refer to these op- 
erators as the cluster operators. Note that we do not 
require that the wave function basis is orthogonal, and 
the nonorthogonality is summarized in an overlap matrix 

^av = Onv{R m — Rn) = {RmH\RnV) ■ 
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FIG. 1. Schematic separation of the elements of a cluster 
dynamical mean field algorithm. 

The next step is to express the Hamiltonian in terms 
of the complete set of operators CR mfJl . In terms of the 
new set of variables it has the form 

H = — t^ v {R m — Rn)CR mtl CR n i> 

+ U ^pd{Rz})c% in c+ 2U c B ^c B , 3P . (1) 

RlnR2VR3PR.it; 

We stress again that the generality of the method. Equa- 



tion (Q) has the form one would obtain by writing the full 
Hamiltonian of electrons in a solid in some tight binding 
non-orthogonal basis. The Hamiltonian is now split into 
three parts, H = H c + H c t, + Ht where H c involves only 
the cluster operators, Hi, contains c# nAI with R n 7^ 
only and plays the role of a " bath" , and finally H c i, con- 
tains both Ctf nU with R and the cluster operators c n . 
Physically H c b couples the cluster with its environment. 
A similar separation can be carried out at the level of 
the action, in the coherent state functional integral for- 
mulation of this problem, where the partition function 
and the correlation functions are represented as averages 
over Grassman variables, 



Z = f II Dc+ R na Dc Rna e 



(2) 
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where the action is given by 
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(3) 



The effective action for the cluster degrees of freedom is 
obtained conceptually by integrating out all the variables 
cr„h with R n ^ in a path integral to obtain an effective 
action for the cluster variables c u , i.e. 



J eff 



-S e/ /[c+ C(1 ] 



/ n D 4, n uDc Rn 



(4) 



R m ^0,n 



Note that the exact knowledge of S e ff allows us to cal- 
culate all the local correlation functions involving cluster 
operators. As described in jjj, this cavity construction 
if carried out exactly would generate terms of arbitrary 
high order in the cluster variables. Our approximation 
neglects the renormalization of the quartic and higher or- 
der terms. Since the action S c b contains only boundary 
terms, the effects of these operators will decrease as the 
size of the cluster increases. Within these assumptions, 
the effective action is parameterized by Go, p „ (r — r'), the 
Weiss function of the cluster and has the form 

Seff = ~ [drdr'J24(T)G^ 1/ (r r')c v (r') (5) 

+ / dT 1 dT 2 dT 3 dT4r t _ lup sC+(T 1 )c+(T 2 )c,;(T4)c p (T 3 ) 
JO 

where r M „ p<; = U uvpi ({0}). Using the effective action 
(||) one can calculate the Green's functions of the cluster 
G c ,^(r-T')[Go\ = -(T t c u (t)c+(t')){Go} and the cluster 
self energies 



Gn 
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(6) 
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b) Self- consistency condition: The cluster algorithm is 
fully defined once a self-consistency condition which in- 
dicates how Go should be obtained from S c and G c is 
defined. In the approach that we propose here the self 
consistent equations become matrix equations expressing 
the Weiss field in terms of the cluster self energy matrix 



A^(iuj n )[e k x, V k \] = ^ - 



kX 



kX,vYkX,v 



(10) 



The hybridization function is related to the Weiss field 
function by expanding Eq. ^ in high frequencies: 



G (iu) n ) = ibj n O - E - A(iun) 



(11) 
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(iu + fi)O(k) - t(k) 



(7) 



with O = 



where 0(k) is the Fourier transform of the overlap ma- 
trix, t(k), is the Fourier transform of the kinetic energy 
term of the Hamiltonian in Eq. (|l|) and k is now a vector 
in the Brillouin zone (reduced by the size of the cluster, 
Lj, in each direction). Equations (||) and (0), can be de- 
rived by scaling the hopping between the supercells as 
the square root of the coordination raised to power of 
the Manhattan distance between the supercells and gen- 
eralizing the cavity construction of the DMFT Q from 
scalar to matrix self energies. If the cluster is defined in 
real space and the self energy matrices could be taken to 
be cyclic in the cluster indices so that the matrix equa- 
tions could be diagonalized in a cluster momentum basis, 
Eq. (^) would reduce to the DCA equation However, 
in the DMFT construction, the clusters have free and not 
periodic boundary conditions, and we treat a more com- 
plicated problem requiring additional matrix inversions. 
c) Connection to the self energy of the lattice: The self- 
consistent solution, G c and E c , of the cluster problem 
can be related to the correlation functions of the orig- 
inal lattice problem through the transformation matrix 
Sn. m a,ia by the equation 



-'lat^aa 



,(k) (8) 



where S is the Fourier transform of the matrix S with re- 
spect to the original lattice indices i. Notice that E; at(J(J / 
is diagonal in momentum and will also be diagonal in the 
variable a if this variable is conserved. 
d) Connection to impurity models: As in single site 
DMFT it is very convenient to view the cluster action 
as arising from a Hamiltonian, 

Himp = ^ t Ep^cj^C^ + ^ ' r fii/p^C^C^ C p C<; 

p? p^p? 
+ ^2 e kxa kX a k x + ^ ( V k^ a k\ c P + h.c.) . (9) 

kX kX,fj, 

Here e k x is the dispersion of the auxiliary band and Vkx,n 
are the hybridization matrix elements describing the ef- 
fect of the medium on the impurity. When the band 
degrees of freedom are integrated out the effect of the 
medium is parameterized by a hybridization function, 



indicating that the impurity 

model has been written in a non-orthogonal local basis 
with an overlap matrix O. 

Let us now consider some examples of this approach. 

a) Single- site DMFT: The simplest example is the single- 
site dynamical mean field theory which is exact in the 
limit of infinite dimensions. In this case the cluster is 
just a single site denoted by 0, and the cluster operators 
are the creation and annihilation operators of that site 
cJ CT , c oa . The cluster Hamiltonian is diagonal in the spin 
variables and reduces to the effective action of the Ander- 
son impurity model. The second step is a scalar equation 
G^ 1 = G~ 1 [G?o]+E c [G ]. Finally the third step identifies 
the self energy of the cluster with the lattice self energy. 

b ) Free cluster: The next example is a free cluster scheme 
for the one band Hubbard model. The method divides 
the lattice into supercells, and views each supercell as a 
complex "site" to which one can apply ordinary DMFT. 
Here R n is the supercell position and a labels the differ- 
ent sites within the unit cell, and the spin. Introducing 
a spin label a and a supercell notation where an atom 
is denoted by the supercell, R n , and the position inside 
the supercell, I; a — (a 1 ',1) and Sn na ^ a — 5 atlT '5ii n +; >ri 
is diagonal in spin and position. In this case the overlap 
matrix is the identity. This real space cluster method has 
been investigated using quantum Monte Carlo methods 
(QMC) by Katsnelson and Lichtenstein Jxcj] . 

c) Multiorbital DMFT in a non-orthogonal basis: An- 
other important special case of our general construction 
is the implementation of single-site DMFT in a non- 
orthogonal basis. In this case the supercell is a single 
site, but the wave functions defining the cluster opera- 
tors are chosen so that they are very localized in real 
space. 

In fact, an implementation of this method, in conjunc- 
tion with a generalization of the interpolative perturba- 
tion theory, as an impurity solver, has resulted in new 
advances in the theory of Plutonium (TTJ Here the flexi- 
bility in the choice of basis is crucial for the success of the 
DMFT program. DMFT neglects from the start interac- 
tions which are not onsite. A high degree of localization 
requires a non-orthogonal basis and the formalism intro- 
duced in this letter. 

d) Other bases: Finally we point out that the most at- 
tractive feature of this method is that it would allow its 
formulation in terms of wave functions which are partially 
localized in real and momentum space such as wavelet 
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functions. This flexibility is most appealing for treat- 
ing problems such as the Mott transition where both the 
particle- like and the wave- like aspect of the electron need 
to be taken into account requiring a simultaneous con- 
sideration of real and momentum space. 

We now prove that the CDMFT approach gives mani- 
festly causal Green's functions. For this we assume, that 
we start the DMFT iteration with a guess for the bath 
function A which is causal. The self energy which is 
generated in the process of solving the " impurity model" 
is also causal. Furthermore, any sensible approximation 
techniques to compute the self energy of the cluster re- 
spects causality, so our proof is valid not just for exact 
solutions of the CDMFT scheme but also for approximate 
solutions as long as the impurity solvers used in the so- 
lution of the cluster impurity problem preserve causality. 
The next step is to show that if a causal self energy is 
introduced in the self-consistency condition, the re- 
sulting bath function A is causal. Since both S c and A 
are matrices, the causality condition needs to be formu- 
lated precisely. For a Fermionic matrix function, A(ui), 
to be causal means that it is analytic in the upper half 
of the complex frequency plane and thus has a spectral 
representation with spectral density =^{A{lu) — A^(u)}, 
and the spectral density matrix is positive definite. It 
is easy to see that the DMFT equations lead to the cor- 
rect analytic properties and the following proof estab- 
lishes the positivity of the bath spectral density. Writing 
£# = 6 — 17 with e, 7 hermitian and 7 positive definite, 
we get: 



(Ajj - A R ) = -2vy 



(12) 



V7 £ 



i + dk 



£ 



-i + d k 



v/7. 



Positivity is reduced to proving that the following matrix 
is negative jl2], i.e. 



2 - 




< 0. 



(13) 



Here wO(fc) - t(k) — e = ^/jdk^/j- Performing a change 
of variables dk — i e ie^J ~it k with 9k a hermitian matrix, 
Eq. ( |i~3| ) reduces to proving that 

K^w-^ + a-^)- 1 ^) (14) 

X \ x \ x ) 

where z = {(e~ 2l9k )) , with (()) denoting an average over 
the Brillouin zone. By performing the substitution |x) = 
(1 - z)\y) Eq. (fil) reduces to 



1 < min 



2 + (y\z + zi\y) 



which clearly holds due to the fact that z is an average of 
unitary matrices, and has the property ||.z^.z|| < 1. From 
Eq. ( [L5| ) and Eq. (^J) it follows that the imaginary part 
of the retarded self energy is always less or equal to zero, 
completing the proof of causality. This proof generalizes 
Rcf. S from scalar to matrix CPA equations. 

In conclusion DMFT has produced a wealth of infor- 
mation in problems where the physics is local and clus- 
ter methods promise to be equally fruitful in more com- 
plex problems where correlations between more sites and 
orbitals need to be taken into account. All the tech- 
niques which have been used for the solution of the sin- 
gle site DMFT are applicable to this cluster extension. 
The most powerful methods of solution have been renor- 
malization group related techniques such as the projec- 
tive self-consistent method Jl3| or the numerical Wilson 
renormalization group techniques [fl4j| . These methods 
carry out a division of the impurity and the bath into a 
low-energy and a high-energy part, and perform an elim- 
ination of the high-energy degrees of freedom in both. 
Extensions of these methods to an abstract cluster is a 
first step in constructing non-perturbative renormaliza- 
tion group method for correlated fermion systems and 
deserves further investigations. 



112/11=1 l + {y\zh\y) + (y\z + z^\y) 
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